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ABSTRACT 

We present the first determination of the intrinsic triaxial shapes and tree-dimensional physical pa- 
rameters of both dark matter (DM) and intra-cluster medium (iCM) for th e galaxy cluster Abell 1689. 
We exploit the novel method we recently introduced (Mo randi et al.|[2010f ) in order to infer the tree- 
dimensional physical properties in triaxial galaxy clusters by combining jointly X-ray and strong 
lensing data. We find that Abell 1689 can be modeled as a triaxial galaxy cluster with DM halo axial 
ratios 1.24 ± 0.13 and 2.37 ± 0.11 on the plane of the sky and along the line of sight, respectively. We 
show that accounting for the three-dimensional geometry allows to solve the discrepancy between the 
mass determined from X-ray and strong gravitational lensing observations. We also determined the 
inner slope of the DM density profile a: we measure a — 0.90 ± 0.05 by accounting explicitly for the 
3D structure for this cluster, a value which is close to the cold dark matter (CDM) predictions, while 
the standard spherical modeling leads to the biased value a = 1.16 ± 0.04. Our findings dispel the 
potential inconsistencies arisen in the literature between the predictions of the CDM scenario and the 
observations, providing further evidences that support the CDM scenario. 

Subject headings: cosmology: observations galaxies: clusters: general galaxies: clusters: individual 
(Abell 1689) gravitational lensing: strong gravitational lensing: weak X-rays: 
galaxies: clusters 



1. INTRODUCTION 

Abell 1689 is a massive cluster with a redshift of 
z = 0.183 and with the largest Einste in radius, around 
45 arcsec for z s = 1, observ ed to date (jBroadhurst et al.l 
l2005t iLimousin et al.l l2007f ) . It has been proposed as a 
standard example of relaxe d object in hydros tatic equi- 
librium (|Xue fc Wul l200l ILemze et al.lT2010i ). but the 
mass derived from the X-ray measurement is twice as 
small as that found from strong gravitational lensing 
at most radii (|Andersson fc Madeiskil 120041: ILemze et al.l 
l2008t iRiemer-Sorensen et all 120091: iPeng et al.l l2009ft. 

iRiemer-Sprensen et al.l (|2009f) showed that the discrep- 
ancy between X-ray and weak lensing mass determina- 
tions has been reduced if we exclude a cool clump plus 
some substructure in the North Eastern part of the clus- 
ter, nevertheless a departure of factor of ~ 2 between SL 
and X-ray mass still remains in the strong lensing region. 

Since lensing is sensitive to the integrated mass con- 
trast along the line of sight, departures from the spheri- 
cal assumption could justify the disagreement between 
X-ray and lensing mass profiles found in the litera- 
ture dGavazzil 2005f) . In particular our recent work 



(|Morandi et al. 120101) presented the first determination 
of the intrinsic shapes and the physical parameters 
of both DM and ICM in the triaxial galaxy cluster 
MACS J1423. 8+2404 by combining X-ray and lensing 
data, and we showed that triaxiality allows to solve the 
long-standing discrepancy between galaxy cluster masses 
determined from X-ray and gravitational lensing obser- 
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vations. Since cluster mass measurements are sensitive 
to the assumptions about symmetry, this suggests that 
clusters with prominent strong lensing features are not 
spherically symmetric and preferentially elongated along 
the line of sight increasing the magnitude of the lensing 
(|Oguri fc Blandfordl[2009t iMeneghetti et al.ll2010D . 

The main motivation for studying Abell 1689 in large 
detail with the best available X-ray and lensing data is to 
get insight into the mass discrepancy and to understand 
if it can be ascribed to standard spherical geometry. In 
this perspective we aim at unveiling its triaxial physical 
properties, shape and mass distribution via implemen- 
tation of the novel method described in iMorandi et al.l 
(|2010[) and by combining X-ray and SL data. 

Hereafter we have assumed a flat ACDM cosmology, 
with matter density parameter Q m = 0.3, cosmological 
constant density parameter J7a = 0.7, and Hubble con- 
stant Ho = 70km/s/Mpc. Unless otherwise stated, we 
estimated the errors at the 68.3 per cent confidence level. 

2. THE DATASET AND THE ANALYSIS 

Here we briefly summarize the most relevant aspects 
of our data reduction and analysis of Abell 1689. 

2.1. Strong lensing analysis 

For the strong lensing analysis we refer to the findings 
of lLimousin et al.l (|2007f ) , who presented a reconstruction 
of the mass distribution of the galaxy cluster Abell 1689 
using detected strong lensing features from deep ACS 
observations and extensive ground based spectroscopy. 
They presented a parametric strong lensing mass recon- 
struction using 34 multiply imaged systems, and they 
inferred two large-scale dark matter clumps, one associ- 
ated with the center of the cluster and the other with a 
north-eastern substructure. We masked out the north- 
eastern sector of both the 2D projected mass map and the 
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X-ray data, in order to avoid the contribution from this 
secondary substructure. We masked out also the central 
30 kpc, which is affected by the mass distribution of the 
cD galaxy. From the lensing analysis the major clump 
of the cluster looks elongated with a minor-major axial 
ratio on the plane of the sky of 1.24 ± 0.13 and position 
angle of 0.4 ± 1 degrees. The 2D projected mass map has 
been rebinned into elliptical annuli, whose eccentricity, 
cent roid and position angle is the same as that inferred 
from lLimousin et al.1 (|2007t ). Then we calculated average 
values of the elliptical symmetric projected mass profile 
k(R), R being the minor radius of the 2D elliptical an- 
nuli, once we masked out the central 30 kpc, which is 
affected by the mass distribution of the cD galaxy. We 
also calculated the covariance matrix C' among all the 
measurements of k(R). 

2.2. X-ray data reduction 

For the X-ray analysis we take advantage of the 2 
Chandra X-ray observations (observation ID 6930 and 
7289) from the NASA HEASARC archive with a total 
exposure time of approx. 150 ks. We summarize here 
the most relevant aspects of the X-ray data reduction 
procedure for Abell 1689. The observations have been 
carried out using ACIS-I CCD imaging spectrometer. 
We reduced these observations using the CIAO software 
(version 4.1.2) distributed by the Chandra X-ray Ob- 
servatory Center, by considering the gain file provided 
within CALDB (version 4.1.3) for the data in VFAINT 
mode. Then we have filtered the data to include the 
standard events grades 0, 2, 3, 4 and 6 only, and there- 
fore we have filtered for the Good Time Intervals (GTIs) 
supplied, which are contained in the fltl.fits file. 
We checked for unusual background rates through the 
lc_sigma_clip, so we removed those points falling out- 
side ±3cr from the mean value. Finally, we filtered ACIS 
event files on energy selecting the range 300-9500 keV 
and on CCDs, so as to obtain an events 2 file. 

2.3. X-ray spatial and spectral analysis 

We outline the methodology of spatial and spectral 
analysis in triaxial galaxy clusters. The general idea is 
to measure the gas density profile in an non-parametric 
way from the surface brightness recovered by a spatial 
analysis, and to infer the observed projected temperature 
profile by a spectral analysis. 

The images have been extracted from the events 2 files 
in the energy range (0.5 — 5.0 keV), corrected by the 
exposure map to remove the vignetting effects, by mask- 
ing out the point sources. We constructed a set of n 
(n = 57) elliptical annuli of minor radius r m around the 
centroid of the surface brightness and with eccentricity 
£;/ (r) fixed to that predicted from the eccentricity eu (r) 
of the DM halo on the plane of the sky from SL data (see 
Sect. 12.41) . The minor radius of each annulus has been 
selected out to a maximum distance i? S pat = 1139 kpc, 
selecting the minor radii according to the following crite- 
ria: the number of net counts of photons from the source 
in the (0.5-5.0 keV) band is at least 200-1000 per annulus 
and the signal-to-noise ratio is always larger than 2. The 
background counts have been estimated from regions of 
the same exposure, which are free from source emissions. 

The spectral analysis has been performed by extract- 
ing the source spectra from n* (n* = 9) elliptical annuli 



of minor radius r* n around the centroid of the surface 
brightness and with eccentricity equal to that predicted 
from the eccentricity (r) of the DM halo from SL data 
(see above). We have selected the minor radius of each 
annulus out to a maximum distance R spC c = 1089 kpc, 
according to the following criteria: the number of net 
counts of photons from the source in the band used for 
the spectral analysis is at least 2000 per annulus and cor- 
responds to a fraction of the total counts always larger 
than 30 per cent. 

All the point sources have been masked out by both 
visual inspection and the tool celldetect, which pro- 
vide candidate point sources. Then we have calculated 
the redistribution matrix files (RMF) and the ancillary 
response files (ARF) for each annulus. 

For each of the n* annul i the spectra have been ana- 
lyzed by using the XSPEC ()Arnaudlll996l version 11.3.2) 
package, by simultaneously fitting absorbed MEKAL 
models m ultiplied by a positive absorption edge as de- 
scribed in lVikhlinin et all ()2005|) to the two observations. 
The fit is performed in the energy range 0.6-7 keV (0.9-7 
keV for the outermost annulus only) by fixing the red- 
shift the redshift at z = 0.183, and the photoelectric 
absorption at the galactic value. We consider three free 
parameters in the spectral analysis for the ith annulus: 
the normalization of the thermal spectrum K[ oc J n\ dV, 
the emission- weighted temperature T* to ^ ; ; the metallic- 
ity Z{ re trieved by employing the solar abundance ra- 
tios from [GrevessTESauSJ (jl998|). Background spectra 
have been extracted from regions of the same exposure, 
which are free from source emissions. 

At last we recover the electron density n e = n e (r] e c >) 
both by deprojecting the surface brightness profile and 
the spatially resolved spectral analysis obtaining a few 
tens (n — 57) of radial measurements in ellipsoidal shells. 
Notice the dependency of n e (r; e c ') on the eccentricity e c ' 
of the ICM along the line of sight to be determined (f or 
further details see Appendix A of lMorandi et al.ll2010f ) . 

The global (cooling-core corrected) temperature T cw 
has been estimated to be T cw = 8.64+q ; 12 keV and 
an abundance of 0.41 ± 0.03 solar value. We classify 
this cluster as a interme diate cooling core source (ICC) 
(jMorandi fc Ettorill2007h : we estimated a t coo \ ~ 3 x 10 
yr. The temperature profile is very regular once we 
masked out the north-eastern quadrant, suggesting a re- 
laxed dynamical state (see upper panel of Fig. [I]). 

2.4. Joint X-ray -hlensing analysis: measuring the 
triaxial physical properties of ICM and DM 

Here we br i efly s ummarize the major findings of 
iMorandi et al.l (|2010( ) for the joint X-ray+lensing analy- 
sis in order to infer tri axial physical properties: fo r fur- 
ther details we refer to IMorandi et al.1 ()2007L I2010Q . 

The lensing and the X-ray emission both depends on 
the properties of the DM gravitational potential well, the 
former being a direct probe of the 2D mass profile and the 
latter an indirect proxy of the 3D mass profile through 
the hydrostatic equilibrium equation applied on the gas 
temperature and density. In order to infer the model pa- 
rameters of both the ICM and of the underlying DM den- 
sity profile, we perform a joint analysis for strong lensing 
and X-ray data. In this perspective we briefly outline 
the methodology in order to infer physical properties in 
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triaxial galaxy clusters: 1) we start with a generalized 
Navarro, Frenk & White (gNFW) triax ial model of the 
DM as described in lJing fe Sutol (j2002D . which is repre- 
sentative of the total underlying mass distribution and 
depends on a few parameters to be determined, namely 
the concentration parameter c, the scale r adius r s and 
the inner sl ope of the DM a 2) following iLee fc Sutol 
(2003l [200l . we recover the gravitational potential and 
surface mass profile fc of a dark halo with such triaxial 
density profile 3) we solve the hydrostatic equilibrium 
equation for the density of the ICM sitting in the gravi- 
tational potential well previously calculated, in order to 
infer a theoretical three-dimensional temperature profile 
T gas in a non-parametric way 4) the joint comparison of 
T gas with the observed temperature and of fc with the ob- 
served surface mass give us the parameters of the triaxial 
DM density model, and therefore all the desired physical 
properties of ICM and DM triaxial ellipsoids (see Fig. 
1). , , 

The work of ILee fc Sutol (120031 ) showed that the ICM 
and DM halos are well approximated by a sequence of 
concentric triaxial distributions with different eccentric- 
ity ratio. We define (eb>) and e c > (e c <) as the eccentric- 
ity of DM (ICM) on the plane of the sky and along the 
line of sight, respectively. The iso-potential surfaces of 
the triaxial dark halo coincide also with the iso-density 
(pressure, temperature) surfaces of the intra-cluster gas. 
Notice that = (ey, u, a) and e c > = e c '(e c >,u,a), 
with u = r/r s , unlike the constant et>', e c i for the adopted 
dark matter halo profile. In the whole range of u, e<j/ e a 
is less than unity (t a je a ~ 0.7 at the center), i.e. the 
intra-cluster gas is altogether more spherical than the 
underlying DM halo. 

In order to infer the model parameters, we construct 
the likelihood performing a joint analysis for SL and X- 
ray data, to constrain the properties of the model param- 
eters q of both the ICM and of the underlying gNFW 
triaxial model of the DM, with 



q = (c,r s ,a,e c 



representing the concentration parameter, scale radius, 
inner slope of the DM and eccentricity of the DM along 
the line sight, respectively. 

The method works by constructing a joint X- 
ray+lensing likelihood: 



C = C X -C 



lens 



(2) 



being £ x and £i ons the likelihoods coming from the X-ray 
and SL data, respectively. 
For £ x holds £ x oc exp(— x x /2), with x x equal to: 



n 

i=l 



(T proj>i (q) - T 



proj,i/ 



(3) 



being T* ro j ; the observed projected temperature profile 
in the ith ring and T UTO \,i ( q) th e convenient projection 
(following Ma zzotta et al.l [2004) of the theoretical 3D 
temperature T(q) recovered by solving the hydrostatic 
equilibrium equation, once we inferred the gas density 
n e (r;e c >) from the brightness and we assume a gNFW 
parametrization for the DM pdm = Pdm(i", q). £icns 



reads: 

_ expj-iKfc^-fc^C-^fc^-fc*)]} 

£lCI1S = (2^)™*/2|C| 1 /2 ' (4) 

where C is the covariance matrix referred to the pro- 
jected mass profile from lensing data including system- 
atic effects (see below), | C| indicates the determinant of 
C, fc* = (fc*, fc^, fc„») are the observed measurements 
of the projected mass profile in the m* elliptical annuli, 
fc(q) the theoretical projected mass profile retrieved by 
our triaxial DM model. 

For the covariance matrix C, the below expression 
holds: 



c 



°sys X 



(5) 



being C' the covariance matrix among the lensing mea- 
surements, X is the identity matrix and cr sys a bias param- 
eter estimator arising from measurements of systematics 
involved in the SL analysis. In order to calculate <7 sys we 
assumed that systematic errors can be described as gaus- 
sian errors via a diagonal matrix cr„„ X with the same 
value in each of the diagonal elements. We checked that 
this simplified assumption docs not affect significantly 
the average value of the physical parameters, while it 
slightly increases (10-20%) their errors. 

We marginalized over (q, <7 sys ) and therefore we have 
£ = £(q,cr sys ). 

So we can determine the physical parameter of the 
cluster, for example the 3D temperature Tj = Tj(q) in 
the jth shell and the elongation e c >(e c r) of the ICM(DM) 
along the line of sight, just by relying on the hydrostatic 
equilibrium equation and on robust results of the hydro- 
dynamical simulations of the DM profiles. In Fig. Q] we 
present an example of a joint analysis for Abell 1689: 
notice that in the joint analysis both X-ray and lens- 
ing data are very well fitted by our model, with a total 
Xtot = Xx + xLs = 7A ( n degrees of freedom), xl = 5.5 
(5 degrees of freedom) and Xiens = 1-9 (2 degrees of free- 



(1) dom), with X\ e 



oc 
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lens ) 



3. RESULTS 



In table Q] we present the best model fit parameters for 
our analysis of Abell 1689. In Fig. [2]we present an image 
of the core of Abell 1689 from optical (Hubble Space Tele- 
scope) observations, with overplotted the projected total 
mass contours computed from the gravitational lensing 
analysis (blue line) and from the X-ray surface bright- 
ness (green line). Our findings outline a picture where 
Abell 1689 is a triaxial galaxy cluster with DM halo axial 
ratios ri' DM b = 1.24±0.13 and r]' BM c = 2.37±0.11, where 
Vbm b i s the axial ratio of the DM on the plane of the sky 
inferred from lensing measurements, and r)' DM c the axial 
ratio of the DM along the line of sight inferred through 
our joint analysis (see table [1]). Notice these elongations 
are statistically significant, i.e. it is possible to disprove 
the spherical geometry assumption. 

The axial ratio of the gas is r? ga s,b' ~ 1.1 — 1-06 (on 
the plane of the sky) and 77 gas .c' ~ 1.6 — 1.3 (along the 
line of sight), moving from the center toward the X-ray 
boundary. 

Here we focus on the implications of our method on 
the CDM scenario and on the discrepancy between X- 
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Fig. 1. — Joint X-ray and lensing analysis for Abell 1689. In the 
upper panel we display the two quantities which enter in the X- 
ray analysis spectral deprojection analysis (eqn. [3j: the observed 
spectral projected temperature Tp ro j ; (big points with errorbars) 
and the theoretical projected temperature T pro j ^(q) (diamonds). 
We also show the theoretical 3D temperature T,(q) (points), which 
generates T pro j i(q) through convenient projection techniques. In 
the lower panel we display the two quantities which enter in the 
lensing analysis (eqn. Q): the observed surface mass profile k* 
(points with errorbars) and the theoretical one fc(q) (solid line). 
The distance of the points for both Tj (q) and k* is representative 
of the true spatial resolution of X-ray and lensing data, respectively. 

ray and lensing masses on Abell 1689, showing that this 
is dispelled if we account explicitly for a 3D geometry. 

3.1. Probing the CDM scenario 

Measuring the three-dimensional mass distribution on 
galaxy cluster scales is a crucial test of the CDM scenario 
for structure formation models, providing constraints on 
the nature of dark matter. Recent works investigat- 
ing mass distributions of individual galaxy clusters (e.g. 
Abell 1689) based on gravitational lensing analysis have 
shown potential inconsistencies between the predictions 
of the CDM scenario relating halo mass to concentration 
parameter, and the relationships as measured in massive 
clusters. 

For example iBroadhurst Sz Barkanal (120081) usi ng the 
distribution of halo profiles from lNeto et al.l (|2007l) found 
that the predicted Einstein radii under the assumption of 
spherical geometry are a factor of two below the observed 
Einstein radii of four massive clusters with spectacular 
Einstein rings (among them Abell 1689). Relatively high 
concentrations parameters of ~ 8 — 14 are derived from 
lensing analysis of A bell 1689 (jBroadhurst et al.l 120051 : 
iLimousin et alll2~007D . These values are larger than the 
concentration parameter expected ba sed on simulation s 
of the standard CDM model (c ~ 4, iNeto et all 120071 ). 
Given that the predicted mass profile is too shallow com- 
pared to the observed ones, the question arises whether 
the projected critical surface density for lensing can be 
exceeded within a substantial radius for this model. 

However, we emphasize the previous analyses em- 
ploy standard spherical modeling of the DM halo, while 



numerical simulations predicts that DM halos show 
axis ratios typically of the order of ~ 0.7 (Sha w et"al"1 
20061). disprov i ng th e spherical geometry assumption. 
Morandi ct al. (2010) demonstrated that the halo triax- 
iality could cause a significant bias in estimating the de- 
sired physical parameters, i.e. concentration parameter 
c, inner slope of the DM a and total mass if a spheri- 
cal halo model is a priori assumed for the model fitting. 
As a consequence, the projected mass distributions of 
the clusters have larger concentration parameter and in- 
ner slope of the DM compared with typical clusters with 
similar redshifts and masses. 

In light of the previous considerations, we evaluated 
the Einstein radius for Abell 1689 via our triaxial joint 
X-ray+SL analysis. The Einstein radius #e occurs at a 
projected radius where the mean enclosed convergence is 
equal to 1. Our triaxial joint X-ray+SL analysis predicts 
9e = 42.7 ± 3.1 arcsec for z s = 1, whi ch is in agreement 
with t he observed value of 45 arcsec (IBroadhurst et al.l 
[20051: ILimousin etldl 120071 iRichard et all 120091) . We 
conclude that the large Einstein radius observed in 
Abell 1689 is not in conflict with CDM predictions, as 
long as the triaxiality of the DM halos are taken into ac- 
count. In this perspective, we also find that the minor- 
major principle axis ratio i]' I)M c = 2.37 ± 0.11 is consis- 
tent wit h the results from numerical simulations within 
- 2.5cr (|Shaw et al.l 120061) . 

Then we focus on the determination of the other pa- 
rameters of the DM halos, namely the concentration pa- 
rameter and the inner slope of the DM. In this per- 
spective, one of the main result of the presented work 
is to measure a central slope of the DM a = 0.90 ± 
0.05 by accounting explicitly for the 3D structure for 
A bell 1689. This value is close to the CDM predictions 
of lNavarro et al.l (jl997l) ( i.e. a — 1), but it is even in bet- 
te r agreement w i th the more recent numeral simulations 
of IMerritt et al.l (|2006l ). which predicts slightly shallower 
inner slope. The value of the concentration parameter is 
4. 58 ±0.34, in agreement with the the oretical expectatio n 
from hydrodynamical simulations of INeto et all (2007), 
where c ~ 4 at the redshift and for the virial mass of 
Abell 1689, and with an intrinsic scatter of ~ 20 per cent. 
This lends support to our insights about the role of the 
effects of geometry on the physical properties and allows 
to solve the arisen potential inconsistencies between the 
predictions of the CDM scenario and the measurements 
in massive clusters. 

If we carry out a standard spherical modeling, we 
obtain the biased value a = 1.16 ± 0.04 for an X- 
ray-only analysis, value larger than that in table [TJ 
The different value of a in triaxial and spherical case 
shows that the systematics involved in neglecting elon- 
gation/flattening of the sources along the line of sight 
are relevant: this likel y justifies the large scatter of a 
found in th e litera t ure (lEttori et alj|2002t IGavazzill2005l: 
Sand et all 120081: IBradac et all 120081: ILimousin et all 
2008: iBiviano fc Saluccill2006l) . 

3.2. Resolving the discrepancy between X-ray and strong 
lensing masses 

Here we briefly summarize the major findings in the lit- 
erature for Abell 1689 in order to study the discrepancy 
between the mass determined from X-ray and gravita- 
tional lensing observations. 
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TABLE 1 

Best model fit parameters of Abell 1689. The columns 1 — 5 refer to the best fit parameters c, r s , a, r\nu,c AND o" sys , while 

THE LAST THREE COLUMNS REFER TO THE MASS AND RADIUS AT A = 2500, RESPECTIVELY, AND TO THE ElNSTEIN RADIUS 



c r s a Vdm,c °"sys M250O -R2500 9 E (z 3 = 1) 
(kpc) (gr/cm 2 ) (10 14 M Q ) (kpc) (arcscc) 



4.58 ±0.34 445 ±35 0.90 ± 0.05 2.37 ±0.11 0.004 ± 0.002 8.58 ± 0.23 556 ± 12 42.7 ±3.1 




Fig. 2. — Optical image of the core of Abell 1689 from the 
NASA/ESA Hubble Space Telescope, with overplotted the pro- 
jected total mass contours computed from the gravitational lensing 
analysis (blue line) and from the X-ray surface brightness (green 
line). The north-eastern sector has been masked out in our joint 
SL+X-ray analysis. 



A recent joint Chandra, HST/ACS, and Sub- 
aru/ Suprime cam analysis by iLemze et al.l (|2008l ) sug- 
gested that the temperature of A1689 could be as high 
as T = 18 keV at 100 hr 1 kpc, almost twice as large as 
the observed value at that radius. The derived 3D tem- 
perature profile was based on the X-ray surface bright- 
ness, the lensing shear, and the assumption of hydro- 
static equilibrium. From the disagreement between the 
observed X-ray temperature and the deduced one, they 
concluded that denser, colder, and more luminous small- 
scale str uctures could bias the X-ray temperature. Nev- 
ertheless iPeng et all (|2009| ) proved that if the temper- 
at ure profile of the a mbient cluster gas is in fact that 
of ILemze et al.l (|2008t ). the cool clumps would have to 
occupy 70-90 per cent of the space within 250 kpc ra- 
dius, assuming that the two temperature phases are in 
pressure eq uilibrium. They con clude that the scenario 
proposed by ILemze et al.l ()2008| ) is unlikely. 

Since lensing is sensitive to the integrated mass con- 
trast along the line of sight, either fortuitous alignments 
with mass concentrations which are not physically re- 
lated to the galaxy cluster or departures of the DM halo 
from spherical symmetry can justify the discrepancy in 
the literature between cluster masses determined from 
X-ray and strong gravitational lensing observat ions, the 
latter being significantly higher than the former (|Gavazzil 
|2005D. 

lLokas et "all (|2006l ) pointed out that Abell 1689 has a 
complex structure in velocity space, suggesting the pres- 



ence of dynamically independent structures along the line 
of sight, whi c h wou ld affect lensing mass estimates, but 
ILemze et al.l (|2009l ) disagree with this projection view 
and they argued that there is no evidence for such sub- 
structure in their velocity data. They conclude that only 
one identifiable substructure at +3000 km/s, 1.50 arcmin 
to the NE, which is seen in the strong lensing mass anal- 
ysis but is determined not to be massive (less than 10 
per cent of the total mass in the strong lensing region). 
Nonetheless, the higher than usual velocity dispersion 
in the cluster center, ~ 2100 km/s, indica tes that the 
central part is quite complex (jCzoskc 2004]). This may 
also imply that the halo is elongated in the line-of-sight 
direction, as galaxies move faster along the major axis. 

When it comes to a "superlens" clusters as Abell 1689, 
halo sphericity is never a justified assumption. Indeed 
lOguri fc Blandfordl (|2009l ) showed that SL clusters with 
the largest Einstein radii constitute a highly biased pop- 
ulation with major axes preferentially aligned with the 
line of sight incre asing the magnitude of lensing, and 
lOguri et all (120051 ) concluded that Abell 1689 weak lens- 
ing measurements are indeed compatible with the CDM- 
based triaxial halo model if Abell 1689 represents a rare 
population (~6 per c ent fraction) of cluster-scale halos. 
iMorandi et al.l (|2010f) demonstrated that triaxiality allow 
to solve the mass discrepancy between the lensing and X- 
ray estimates in the galaxy cluster MACS J1423. 8+2404. 

Indeed, in Fig. [3] we compare the 2D mass enclosed 
within circular apertures of radius R for lensing, for an 
X-ray-only analysis under the assumption of spherical 
geometry, and for a joint X-ray+lensing analysis tak- 
ing into account the 3D geometry. We emphasize the 
good agreement between the masses inferred from lens- 
ing and a joint analysis based on triaxial modeling. On 
the contrary an X-ray-only analysis based on the stan- 
dard spherical modeling clearly predicts systematically 
lower masses by ~ a factor of two in the radial range 
between 30 and 400 kpc. 

This confirms our insights about the role of the effects 
of geometry on the physical properties and solve the long- 
standing discrepancy between X-ray and strong lensing 
mass of Abell 1689. 

4. SUMMARY AND CONCLUSIONS 

In this paper we have employed a triaxial halo model 
for the galaxy cluster Abell 1689 to extract more reliable 
information on the three-dimensional shape and physi- 
cal parameters, by combining X-ray and strong lensing 
measurements. 

We demonstrated that the halo triaxiality could cause 
a significant bias in estimating the desired physical pa- 
rameters, i.e. concentration parameter c, inner slope of 
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R (kpc) 

Fig. 3. — 2D masses enclosed within a circular aperture of radius 
_R from lensing data (points with errorbars), from an X-ray-only 
analysis under the assumption of spherical geometry (solid line with 
the 1-(T error gray shaded region), and from a joint X-ray+lensing 
analysis taking into account the 3D geometry (dot-dashed line with 
the l-cr error gray shaded region). 
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the DM a and total mass if a spherical halo model is a 
priori assumed for the model fitting. 

We focused on the implications of our method on the 
CDM scenario, proving that the value of the c and a are 
in agreement with the CDM predictions, once we prop- 
erly accounted for the 3D shape of the cluster. Depar- 
tures of c and a from the theoretical expectation of the 
CDM scenario found in the literature can be explained 
by a halos having the major axis preferentially oriented 
toward the line of sight. In particular, accounting for 
the 3D geometry allows to resolve the long-standing 
discrepancy between X-ray and strong lensing mass of 
Abell 1689 in literature and predicts an Einstein radius 
in agreement with the observations. 
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